Shot noise and Coulomb blockade of Andreev reflection 
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We derive low energy effective action for a short coherent conductor between normal (N) and 
superconducting (S) reservoirs. We evaluate interaction correction 8G to Andreev conductance and 
demonstrate a close relation between Coulomb effects and shot noise in NS systems. In the diffusive 
limit doubling of both shot noise power and charge of the carriers yields \8G\ four times bigger than 
in the normal case. Our predictions can be directly tested in future experiments. 



It is well known that low energy electron transport 
across the interface between normal metals and super- 
conductors (NS) is provided by the mechanism of An- 
dreev reflection 1 . This mechanism involves conversion of 
a subgap quasiparticle entering the superconductor from 
the normal metal into a Cooper pair together with simul- 
taneous creation of a hole that goes back into the normal 
metal. Each such act of electron-hole reflection corre- 
sponds to transferring twice the electron charge e* = 2e 
across the NS interface and results, e.g., in non-zero con- 
ductance of the system at subgap energies 2 . 

Let us assume that two bulk metallic electrodes, one 
normal and one superconducting, are connected by an 
arbitrary - though sufficiently short - coherent conductor 
as it is schematically shown in Fig. 1. This conductor is 
characterized by the normal state conductance 



N 



mm 



FIG. f : Short coherent conductor connecting normal and su- 
perconducting reservoirs. 



well known form 7 



(\6I\ 2 )=e\V\G N N , Pn = 
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where T n define transmissions of all conducting channels 
and the factor 2 accounts for spin degeneracy. Evaluat- 
ing the conductance Ga of the NS structure in Fig. 1, 
at temperatures/voltages well below the superconducting 
gap A one finds 2 - 



where V is the average voltage across the conductor. In 
NS systems Andreev reflection also leads to the current 
shot noise at energies below the superconducting gap. In 
this case in the zero energy/frequency limit and at T — > 
one obtains^ 



<|<5/| 2 ) = 2e\V\G A p A , f3 A = 



E n T n (l-T n ) _ 

Tn Tn 
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where "Andreev transmissions" T n are related to T n as 



% = Til '(2 - T n f 



(3) 



Comparing Eqs. (J5J), ([3]) with the Landauer formula (JTJ 
one immediately observes that Andreev conductance G a 
formally describes "normal" transport of spinless quasi- 
particles (hence, no extra factor 2 in front of the sum) 
with charge e* = 2e across some effective coherent scat- 
terer with channel transmissions T n ([3|). 

Later it was realized that this formal analogy applies 
not only to electron transport but also to low frequency 
shot noise^^ and eventually to full counting statis- 
tics (FCS)&. Consider, for instance, current fluctuations 
51 (t) = I(t) — I around its average value I = (l(t)). 
In normal conductors at T — > and in the zero fre- 
quency limit the correlator for such fluctuations has the 



where T n is again defined by Eq. ([3]). Again, a close 
similarity between Eqs. ((H) and |(SJ) is obvious: The re- 
sult iJSJ) just describes shot noise produced by carriers 
with effective charge e* = 2e in a coherent scatterer with 
conductance Ga and Fano factor Pa- 
Id. the important case of diffusive NS structures dou- 
bling of the carrier charge also implies doubling of the 
shot noise 4 . In this case the sums over transmission 
channels in the above equations can be evaluated in a 
straightforward manner with the results 



Gn = Ga, Pn — Pa = 1/3, 



(6) 



which yield (\SI\ 2 ) = 2e\V\G N /3 for NS structures and 
the two times smaller result in the normal case. This 
doubling of the shot noise in diffusive NS systems was 
indeed observed in experiments^. 

More recently another interesting observation was 
reported^. The authors of this experiment investigated 
short metallic nanowires attached to bulk superconduct- 
ing electrodes. In a number of samples superconductiv- 
ity inside the wire was destroyed due to phase slippage 
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and, hence, such samples effectively represented hybrid takes the form 
normal-superconducting structures, e.g., similar to those 
depicted in Fig. 1. Remarkably, the authors 10 discovered 
that as long as the electrodes stayed superconducting the 
measured I-V curves could be well fitted by the theory 
of Coulomb blockade in normal coherent conductors 11 
provided the electron charge e was substituted by some 
effective charge q c g larger than e but smaller than 2e. If, 
however, superconductivity in bulk electrodes was sup- 
pressed, the I-V curves of exactly the same form but 
with q c g ~ e were observed. Although these observa- 
tions strongly indicate that Andreev reflection can be 
involved, no theoretical explanation of the experiments 10 
was offered until now. 

Below we develop a theory describing an interplay be- 
tween Coulomb blockade and Andreev reflection. We will 
explicitly evaluate the interaction correction to Andreev 
conductance and demonstrate its direct correspondence 
to the shot noise in NS structures. Hence, very different 
experiments^ 9 - andi - turn out to be closely related. Both 
measure the same effective charge, i.e. e* = q c ff. 

Effective action. As it is shown in Fig. I, we will 
consider big normal and superconducting reservoirs con- 
nected by a rather short normal bridge (conductor) with 
arbitrary transmission distribution T n of its conducting 
modes (for each n the value T n is the same for spin-up 
and spin-down electrons). Both phase and energy relax- 
ation of electrons may occur only in the reservoirs and 
not inside the conductor which length is assumed to be 
shorter than dephasing and inelastic relaxation lengths. 
In contrast to^ (where the Thouless energy exh of the 
normal conductor plays an important role), here exh of 
the bridge is irrelevant as it is supposed to be higher 
than any other energy scale in our problem. As usually, 
Coulomb interaction between electrons in the conductor 
area is accounted for by some effective capacitance C. 

In order to analyze electron transport in the presence 
of interactions we will make use of the effective action for- 
malism combined with the scattering matrix technique. 
This approach - very successful in the case of normal 
conductor s 11 ! 13 ' 14 ^ 5 - can be conveniently generalized to 
the superconducting case. In fact, the structure of the ef- 
fective action remains the same also in the latter case, one 
should only replace normal propagators by 2 x 2 matrix 
Green functions which account for superconductivity, as 
it was done, e.g., i n 16 i 17 ' 18 . 

Following the usual procedure we express the kernel J 
of the evolution operator on the Keldysh contour in terms 
of a path integral over the fermionic fields which can be 
integrated out after the standard Hubbard-Stratonovich 
decoupling of the interacting term. Then the kernel J where 



J = J T>tpiT>tp2 exp(iS[tp]), 



(7) 



where <p\ t 2 are fluctuating phases defined on the for- 
ward and backward parts of the Keldysh contour and 
related to fluctuating voltages 1^,2 across the conductor 
as ifi,2(t) — eV\ t 2- Here and below we set h = 1. 

The effective action consists of two terms, S[<p] = 
S c [cp] + S t [<p] , where 



iS c [V] = ^ J dt'&l = % J dt( P + V>~ (8) 



describes charging effects and the term St [V] accounts for 
electron transfer between normal and superconducting 
reservoirs. It readsi& 



^Trln 



T, 



l + ^({G N ,G s }-2) 



(9) 



where Gn and Gs are 4x4 Green-Keldysh matrices of 
normal and superconducting electrodes which product 
implies time convolution and which anticommutator is 
denoted by curly brackets. In Eq. © we also introduced 
"classical" and "quantum" parts of the phase, respec- 
tively ip+ — ((pi + </52)/2 and <p~ = <pi — <p2- 

Without loss of generality we can set the electric po- 
tential (and, hence, fluctuating phases) of the supercon- 
ducting terminal equal to zero. Then the Green-Keldysh 
matrix of this electrode can be written in a simple form 



G R/A (t) = ±5{t)f 3 - 9{±t) (Af 3 Ji(A*) + zAJ (Ai) 



and G K = G R F - FG A , where F(t) = -iT/ sinh[7rTt] 
is the Fourier transform of I — 2n(e) and n(e) = 1/(1 + 
e e / T ) is the Fermi function. Here Jo,i are the Bessel 
functions, f» are the Pauli matrices, 0(f) is the Heaviside 
step function and A = iAf 2 , where A is chosen real. 

The Green-Keldysh matrix of the normal terminal is 
defined as 



& N (t,t') = U\ K ]Q.v(/./') 
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oitpi(t)T 3 







iip2(t)f 3 



(1 - 2n(e))T 3 2n(e)f 3 
2(l-n(e))f 3 (2n(e)-i; 
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-i<P2(t')f 3 



(11) 
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Substituting the above expressions for Gs and Gn into 
Eq. © we arrive at the action which fully describes 
transfer of electrons between N- and S-terminals to all 
orders in T n . 

In the limit of small channel transmissions one can 
expand St in powers of T n . Keeping the terms up to ~ T"„ 
one recovers the contribution from Andreev reflection. At 
low energies this part of the action reduces to the same 
formii as that for normal tunnel barriers^ in which one 
substitutes e by 2e and Gn by G a- Here, however, we 
are aiming at a more general description which includes 
arbitrary transmission values T n . For this reason we will 
proceed differently. 

Let us define the matrix Xo [<£>+] = 1 — T n /2 + 
(T n /A){G N ,G S }\ l p_=o. As the action St vanishes for 
(f-(t) = one has TrlnAo = 0. Making use of this 
property we can identically transform the action @ to 
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Trln [l + Xq 1 o X'] 



(12) 



where X' = l + (T„/4) ({G N ,G S } - 2)-X . At temper- 
atures and voltages well below the superconducting gap 
Andreev contribution to the action dominates. Hence, it 
suffices to consider the limit of low energies e« A and 
set G s ->■ ( n 2 2 ) • Then we obtain 



2-T„ 

S(t, t')l -^h sin [<p + (t) - y+(t')] F{t, t')ih 



(13) 







<J(t,t')l 



and 



X'(t,t') = ^S(t,t') 






— sin </3_ (t)iT2 


smip-(t)ii~2 






-T n F{t,t') 



- cos v ^ sin v £ ^ cos{tp + (t) - tp + (t'))if 2 




COS y 2 ' COS y £ ~ 1 


sin(<p + (i) - tp + (f))iT 2 


sin ' £ ^f L sin ^ sm(tp + (t) - <p + (t'))iT 2 


sin y ^ cos v ^ 1 cos(ip + (t) — Lp + (t'))iT 2 



(14) 



Now let us assume that either dimensionless Andreev 
conductance gA = ^^2 n %i is large, > 1, or tem- 
perature is sufficiently high (though still smaller than 
A). In either case one can describe quantum dynam- 
ics of the phase variable (p within the quasiclassical 
approximatio n 11 ' 13 which amounts to expanding St in 
powers of (small) "quantum" part of the phase tp-(t). 
Employing Eqs. (j 1 2[) - (Tl4|) and expanding St up to terms 



ip_ we arrive at the Andreev effective action 
iS t = iS R - Si, 



(15) 



where 



iSf, 



' 2tt 



dt' (p-(t')cp + (t'), 



Si 



-— I dt 1 I dt" 



(16) 
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4 J ' J ' sinh 2 [7rT(^-i")] 
o o 

x[/3 A cos(2^+(t') - 2if>+{t")) + 1 -p A }. (17) 

Eqs. (fT5]) - (pT7|) represent the central result of our work. 
It is remarkable that the action St is expressed in exactly 
the same form as that for normal conductor s 11 ' 13 derived 
within the the same quasiclassical approximation for the 
phase variable ip(t). In order to observe the correspon- 
dence between the actio n 11 ' 13 and that defined in Eqs. 
(Ti"5]) - (|17[) one only needs to interchange 
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Ga, Pn *-> Pa 



(18) 



and to account for an extra factor 2 in front of the phase 
ip + under cos in Eq. (|17p . This extra factor implies 
doubling of the charge during Andreev reflection. 

Shot noise and interaction correction. Further analy- 
sis is formally similar to that o£ii. Hence, we can im- 
mediately proceed to our final results. Let us define the 
average current and the current-current correlator as 



(/(*)) = ie / V<p± 



iS[<p] 



:(H)- 



■6ip-(t) 
V<p± 



5ip-(t)6<p-(t') 



(19) 
iS &\ (20) 



where (77)+ = {I(t)I{t') +I(t')I(t)). In the absence of 
interactions we set tp+ = eV and trivially recover the 
standard result I = GaV. For the current fluctuations 
SI(t) from Eqs. p3|) - ([2H|) analogously toii we obtain 



(\si\l) 

G A 



(1 



/^Vcoth — 



Pa \-v , „ Trt , wi 2eV 

-^-^( UJ ±2eV) coth 2T . 



(21) 



This equation fully describes current noise in NS struc- 
tures at energies well below the superconducting gap. For 
eV > T, uj Eq. |J2I]) reduces to the result^ (J3J) while in 
the diffusive regime the correlator ([2Tj) - together with 
Eqs. $Q - matches with the semiclassical result^. 

Let us now turn on interactions. In this case one should 
add the charging term to the action and account for 
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FIG. 2: The interaction correction SG — dl / dV —G n for short 
diffusive conductors at T = Gn/^-kC. The upper and lower 
curves correspond to normal and NS structures respectively. 



phase fluctuations. Proceeding along the same lines as 
for g A > 1 or max(T, eV) > E c = e 2 /2C we get 



I = G A V- 2e(3 A Tlm 



1 



ID 



(22) 



where ^(x) is the digamma function, w — qaEc /7t 2 T- 
and v = 2eV/7rT. This result is plotted in Fig. 2. 

The last term in Eq. (|2"2"|) is the interaction correction 
to the TV curve which scales with Andreev Fano factor 
(3a in exactly the same way as the shot noise (0. Thus, 
we arrive at an important conclusion: interaction correc- 
tion to Andreev conductance of NS structures is propor- 
tional to the shot noise power in such structures. This 
fundamental relation between interaction effects and shot 
noise goes along with that established earlier for normal 
conductor s 1 ^ 20 extending it to superconducting systems. 
In both cases this relation is due to discrete nature of the 
charge carriers passing through the conductor. 



Another important observation is that the interaction 
correction to Andreev conductance defined in Eq. (|22|) 
has exactly the same functional form as that for normal 
conductors, cf. Eq. (25) inii. Furthermore, in a spe- 
cial case of diffusive systems due to Eqs. © the only 
difference between the interaction corrections to the I-V 
curve in normal and NS systems is the charge doubling in 
the latter case. As a result, the Coulomb dip on the TV 
curve of a diffusive NS system at any given T is exactly 2 
times narrower than that in the normal case. We believe 
that this narrowing effect was detected in normal wires 
attached to superconducting electrodes^, cf. Fig. 3c in 
that paper—. 

The above discussion demonstrates that seemingly dif- 
ferent experiments^ andi - are actually closely related: 
Doubling of the shot noise found in NS structures^ cor- 
responds to narrowing of the TV curves observed inia, i.e. 
e* = q e g. The key reason behind this correspondence 
is the relation between shot noise and interaction cor- 
rection to conductance in NS systems established above. 
The absolute value of this interaction correction is pro- 
portional to (effective charge) x (shot noise power), i.e. 
doubling of the shot noise in diffusive NS structures im- 
plies 4 times bigger interaction correction to conductance 
than in the normal case, see Fig. 2. The above predic- 
tions can be verified by independently measuring shot 
noise and Coulomb blockade effects in the same NS struc- 
ture, e.g., as it was already done in normal conductors 22 . 

In summary, we theoretically described the interplay 
between Coulomb blockade and Andreev reflection and 
demonstrated a direct relation between shot noise and 
interaction effects in NS systems. Further extension of 
our theory will include the impact of interactions on FCS. 

This work was supported in part by RFBR grant 09- 
02-00886. 
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